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ABSTRACT 

Despite the great variation in the light curves of Gamma Ray Burst (GRB) prompt emission, their 
spectral energy distribution is generally curved and broadly peaked. In particular, their spectral 
evolution is well described by the hardness-intensity correlation during a single pulse decay phase, 
when the SED peak height Sp decreases as its peak energy Ep decreases. We propose an acceleration 
scenario, based on electrostatic acceleration, to interpret the Ep distribution peak at ^ 0.25 MeV. 
We show that during the decay phase of individual pulses in the long GRB light curve, the adiabatic 
expansion losses likely dominate the synchrotron cooling effects. The energy loss as due to adiabatic 
expansion can also be used to describe the spectral evolution observed during their decay phase. The 
spectral evolution predicted by our scenario is consistent with that observed in single pulses of long 
BATSE GRBs. 

Subject headings: stars: gamma-ray burst: general, radiation mechanisms: non-thermal, acceleration 
of particles. 



1. INTRODUCTION 

The acceleration mechanisms and the radiative pro- 
cesses underlying the prompt emission in long GRBs 
are still unclear. Their large energies released on short 
timescales likely requires that the radiation is produced 
in a highly relativistic jet (e.g. Mcszaros 2002). 

We may consider that a GRB consists of fundamen- 
tal units of emission, or pulses, in the light curve (e.g. 
Norris et al. 1996, Stern & Svensson 1996). The pulse 
structures generally show a sharp rise and a slower decay 
phase to the background flux threshold. However, GRBs 
exhibit a wide variety of light curves, both in shape and 
in duration, and the combination of many such pulses 
could create the observed diversity and complexity of 
light curves (Fishman et al. 1994). 

GRB spectral energy distributions (SEDs) are gener- 
ally curved and broadly peaked. They are usually well 
described by the Band function (Band et al. 1993) or, 
as shown more recently, by the simpler (and physically 
motivated) log-parabolic model (Massaro et al. 2010, 
hereinafter MIO). Kaneko et al. (2006, 2008) investi- 
gated the time resolved spectral behavior of the GRBs 
in the BATSE catalog (Paciesas et al. 1999) and showed 
that the distribution of their SED energy peaks, Ep, is 
symmetric around ~ 0.25 MeV (Goldstein et al. 2010). 

As reported in the analysis of Ryde & Svensson (2002), 
the GRB spectral evolution during the decay phase of 
individual pulses can be described by the hardness- 
intensity correlation (HIC) between the time resolved 
SED peak height, Sp (defined as vF^, and proportional to 
the total flux measured at Ep), and the peak energy Ep 
in the form of a power-law: Sp oc E^, where the distribu- 
tion of the 1] parameter peaks at value 1.7 (Borgonovo 
& Ryde 2001). 

Subsequently, Ryde & Petrosian (2002) showed that 
a powcrlaw form (i.e. Sp oc E^) can be reproduced 
through kinematic effects when applied to a spherical 
shell expanding at extreme relativistic velocity. The cur- 



vature of a relativistic shell would make the photons emit- 
ted off the line of sight delayed and affected by a varying 
Doppler boost as a result of the increasing angle at which 
the photons are emitted with respect to the observer. 
They show that these so called curvature effects, charac- 
terized by a time scale Tang, display a similar trend to 
that of the hardness-intensity correlation (HIC) observed 
in the GRB spectral evolution, with the parameter rj =2 
then expected. However, they also argued that an intrin- 
sic correlation between these two spectral parameters in 
the GRB prompt emission could affect the observed HIC 
(see also Kocevski et al. 2003). These curvature effects 
are dependent on the radius R of the emitting shell, and 
are likely to be negligible if _R < 10^"* cm. 

Motivated by these observations, wc propose an accel- 
eration scenario to explain the observed Ep distribution 
of individual pulses in long GRBs. We also show that 
losses for adiabatic expansion play an important role dur- 
ing the decay phase and could be more relevant than syn- 
chrotron radiative cooling. Finally, wc argue that the ob- 
served spectral behavior of long pulses, during the GRB 
decay phase of the light curve, can be described taking 
into account the energy loss for adiabatic expansion. 

For our analysis, we use cgs units and we assume a flat 
cosmology with _ffo = 72 km s~^ Mpc"-'^, VIm = 0.26 and 
= 0.74 (Dunkley et al. 2009). Unless stated oth- 
erwise, primed quantities refer to the observer reference 
frame while unprimed quantities refer to the GRB frame. 

2. ELECTROSTATIC ACCELERATION MECHANISM 

We propose a particle acceleration scenario to explain 
the Ep distribution around the observed value of ~ 0.25 
MeV. We assume that the acceleration mechanisms oc- 
curring during the GRB prompt emission are a combina- 
tion of systematic acceleration, responsible for the energy 
peak position of the accelerated particle energy distribu- 
tion (PED), and stochastic acceleration, which accounts 
for the broadening of the PED around its peak (MIO). 

As proposed by Cavaliere & D'Elia (2002) for Blazar 
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jets, GRBs could be powered by the Blandford & Znajek 
mechanism (1977) for the extraction of rotational energy 
from a spinning Black Hole (BH) via the Poynting flux 
associated with the surrounding magnetosphere. In these 
magnetospheres, the electric fields parallel to magnetic 
fields can accelerate charged particles. They can arise, 
for example, as a result of magnetic field reconnection in 
current sheets or MHD jet instabilities (e.g. Litvinenko 
1996, Medvedev & Loeb 1999). 

The force free condition E x B = governing these 
magnetospheres breaks down when the electric field E < 
B. In particular, electric fields are electrodynamically 
screened out at distances that exceed the Debye length, 
d, that for a pair plasma is defined as: 
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where ujp is the plasma frequency, 7 is the electron 
Lorentz factor, me is the electron mass, e its electric 
charge, c the speed of light, B the magnetic field and n 
the plasma density. Electric fields parallel to magnetic 
fields accelerate charged particles and consequently, the 
particle energy gain for each acceleration step can be 
written as: 
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Substituting d from Equation (1), we obtain an expres- 
sion for the Lorentz factor of the accelerated particle: 
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We note that the above expression is similar to the as- 
sumption that the electron energy density Ue ~ n 7 m 
is twice the magnetic energy density ub = B'^/Sn, close 
to the equipartition condition. 

With the above acceleration scenario, for an emitting 
region with particle density n ~ 5x10^ cm^^, a mag- 
netic field B ~ 10'' G and a beaming factor S ~100, all 
of which are typical values for GRB models (e.g. Zhang 
& Meszaros et al. 2002), the synchrotron energy peak E'^ 
is ~ 0.3 MeV, in agreement with the observed E'p distri- 
bution. We argue that the variance of the E'^ distribution 
can be due to the dispersion of the other parameters and 
their intrinsic variations during the burst. The Poynting 
flux in the current sheet provides its luminosity, which 
can be estimated as L = cB'^/2n x {Ics w), where w and 
Ics are the current sheet width and its length, respectively 
(Litvinenko 1999). 

The typical observed isotropic luminosity of a GRB 
is -^10'''^ erg s^^ so the intrinsic equivalent value. 



rescaled for a beaming factor of 100, is Li. 



~ 10*^ erg s ^. Assuming Z^s ~ w ~ 10^^ cm, as de- 
rived from the GRB variability timescale (i.e. ~ 0.1 s), 
the Poynting flux in a single current sheet is L ~ 10^** 
erg s~^, the same order of magnitude of the GRB intrin- 
sic luminosity. 

3. PARTICLE ENERGY LOSSES 

A simple scenario to describe single pulses in long GRB 
light curves assumes an impulsive heating of particles 
and a subsequent cooling and emission. The rise phase of 
pulses is attributed to particle acceleration energizing the 



emitting region while the decay phase reflects the particle 
energy losses. In the following, we show that adiabatic 
expansion is the main process responsible for the parti- 
cle energy losses during the decay phase of single pulses. 
This is also supported by the fact that the synchrotron 
cooling time appears too short to account for the decay 
of GRB pulses. In addition, the observational evidence 
that GRB SEDs are curved (e.g. log parabolic) and not 
the superposition of two power laws (e.g. Band func- 
tion) is a strong indication that stochastic acceleration 
occurs during the prompt emission (MIO). This suggests 
that both systematic acceleration (e.g. due to electric 
fields) and stochastic acceleration mechanisms (e.g. due 
to turbulence) balance the synchrotron radiative losses. 

We neglect the radiative losses from inverse Compton 
emission, since GRB prompt emission does not appear to 
be dominated by the high energy 7-ray component (i.e. 
> 100 MeV, Abdo et al. 2009). 

The adiabatic expansion of the emitting region can be 
described by a self similar model in which the temporal 
evolution of the radius and consequently the density, can 
be expressed as: 



R{t) = Ro 



, n{t) = no 



-3q 



, (4) 



where q is the expansion index and is positive (i.e. 
< (J < 1), and to is the reference time. The rate of 
expansion of the emitting region is defined by the rela- 
tion: 
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and the adiabatic expansion losses of a single particle can 
be written as: 
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for which the analytical solution is: 
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Assuming conservation of the magnetic flux, we can ex- 
press the temporal evolution of its value as: 
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where Bq is the initial value at time to. We note that 
by replacing the relation for the magnetic flux conserva- 
tion (Equation (8)) and for the temporal evolution of the 
particle density (Equation (4) in Equation (3), we obtain 
the same temporal evolution for the particle energy as 
derived in Equation (7). Consequently, the assumption 
of the conservation of the magnetic flux is in agreement 
with the electrostatic acceleration scenario. 

The ratio p between the energy losses due to syn- 
chrotron radiative cooling and adiabatic expansion is: 
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where Bq is the magnetic fleld evaluated at the initial 
time to, and cft is the Thomson cross section. 



Gamma Ray Burst prompt emission 
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By substituting Equation (3) into Equation (9), the 
ratio of energy losses p can be also written as: 

p= 1.27x10-4 ^ flV (10) 



q no 



We assume q > 1/5, as in the case of constant expansion 
rate (i.e. q = 1) or expansion with constant energy (i.e. 
q = 2/5, Sedov phase condition). Then, for typical values 
of the prompt emission region, during ther decay phase, 
Bq ~ 10^ G and n ~ 10'' cm~^, the synchrotron radiative 
losses are dominated by those of the adiabatic expansion 
(i.e. p < 1) after a time interval of ^ 1 s. We also 
note that for a more compact source with Bq ~ 10'' G 
and n ~ lO^cm is ^ 10 s. 

Finally, we remark that the synchrotron cooling time 
is dependent on the particle energy, with Tsyn — 5.16 x 
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s. In contrast, the energy loss time due to 



adiabatic expansion. Tad = l/iad — R{t)/R (see Equa- 
tion (6), is energy independent and so does not affect the 
shape of the particle energy distribution during the decay 
phase. Therefore, if the synchrotron cooling is the main 
energy loss mechanism, we expect a drastic change in the 
curvature of the observed spectra, making them narrower 
with time, while energy losses for adiabatic expansion 
should preserve the spectral curvature (e.g. Massaro et 
al. 2006). 

4. SPECTRAL EVOLUTION OF PULSE DECAY PHASE 

To interpret the hardness-intensity correlation we note 
that in a synchrotron scenario, the intrinsic energy peak 
of the spectral energy distribution depends on the mag- 
netic field and the electron energy as iJp ~ 7^ B, while 
the maximum of the SED evaluated at this energy is Sp ^ 
N j'^ B^, where N is the total number of emitting parti- 
cles (e.g. Massaro et al. 2008). Applying the corrections 
due to the adiabatic expansion losses (i.e. Eq. 7 and 
Eq. 8) and assuming the number of emitting particles 
is (approximately) constant, we find that the intrinsic 
parameters Ep and 5*^ have temporal dependence: 
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This gives the expected intrinsic relation between the 
two SED parameters Ep and Sp as a power-law: Sp oc 
{EpY^"^, independent of the value of the expansion in- 
dex q. The relativistic corrections (due to the relativis- 
tic beaming or curvature effects, e.g. Ryde & Petrosian 
2002) do not affect the intrinsic correlation between Ep 
and Sp. Therefore, the expected observed power-law in- 
dex is still ~ 1.5, and thus near the peak of the ry distri- 
bution estimated for the pulse decay phase of long GRBs. 
In addition, assuming that the size of the emitting region 
is ^10^'^ cm, the curvature effects are negligible, because 
their timescale is too short Tang 10~^ s to explain the 
pulse decay in long GRBs. Thus, we describe the de- 
cay phase in the GRB prompt emission has energy loss 
dominated by adiabatic expansion, assuming that the ac- 
celeration energy gain, via systematic and stochastic ac- 
celeration, balances the synchrotron radiative cooling.. 

We note that the PED in the form of a log-parabolic 
function is a good approximation for the solution of the 
kinetic equation for the particles when considering terms 



taking into account systematic and stochastic accelera- 
tion as well as including synchrotron radiative losses and 
adiabatic expansion (Kardashev 1962, Tramacere et al. 
2009, Paggi et al. 2009, MIO). 
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Fig. 1. — The spectral evolution (from the thick red line to the 
dashed one) of the decay phase for an adiabatically expanding syn- 
chrotron shell based on the assumption described in Sect. 3. Each 
SED is evaluated every 0.5 s for a decay phase of 6 s. It is evident 
that the shape of the SED is preserved and the curvature does not 
change. In our calculations we set Bq = 10^ G, no = 10^ cm~^, 
7p is 2x10*, 5 = 50 and q = 0.8. The PED spectral curvature is r 
= 5.0 and the initial volume to 10"^^ cm^. 

We calculate the synchrotron emission of an adia- 
batically expanding spherical region assuming an emit- 
ting PED with a log-parabolic shape: TV (7) = 
^0 (7/7p)-^-''^°s(''/'^''^. where 7p (i.e. the peak of 7^ x 
N{'-f)) is (7^)^^^ , r is PED curvature and n (x Nq is the 
density (see MIO). Then, as shown in Figure 1, the spec- 
tral curvature of the synchrotron SED is constant during 
the pulse decay phase. 



++++++ 



Fig. 2. — The spectral curvature, measured during the v^hole 
light curve of the long, single pulses in GRB 941026 (top panel) 
and GRB 950818 (lower panel), is seen to be relatively constant. 
The time resolved spectral analysis has been performed with the 
log-parabolic function (see MIO). The shape of the light curve 
is given by dashed line; right hand flux scale). The background 
BATSE light curve (25 - 55 keV) is also shown (dashed line). 

In Figure 2, we show the time resolved spectral analysis 



4 



F. Massaro & J. E. Grindlay 2010 



of the decay phases for two long, single pulse BATSE 
GRBs, to compare our model with the observed spectral 
behavior. 

We adopted the log-parabolic function (i.e. F{E) — 
Sp /E^ {E/Ep)-''^''3iE/E^)^ I^;^Q) describe the time 
resolved SED and to measure the spectral curvature b. 

In Figure 2, we show that, the examples of GRB 941026 
and GRB 950818 do not show significant variation of 
their curvature b during their pulse decay phases, in 
agreement with the scenario dominated by adiabatic ex- 
pansion losses. 

The detail of the spectral behavior observed during the 
decay phase (i.e. between 66 and 70 s) of GRB 950818 
is also shown in Figure 3. The curvature is not drasti- 
cally varying over the whole burst, in agreement with our 
scenario (Figure 1). 

We also note the presence of small fluctuations in b, 
which appear to be anti-correlated with secondary peaks 
in the GRB light curve. As already pointed out by Vet- 
ere et al. (2007) from the analysis of the BeppoSAX 
WFC archive, the low energy (i.e. 2 - 30 keV) GRB light 
curves are characterized by peaks superposed on a slowly 
evolving component. A more detailed analysis of these 
fluctuations will be presented in a forthcoming paper. 

5. CONCLUSIONS 

We propose an electrostatic acceleration scenario to in- 
terpret the Ep distribution of GRB time resolved SEDs. 
We show that taking into account adiabatic expansion 
losses it is possible to describe the spectral evolution dur- 
ing the decay phase of individual pulses in long GRBs. 
Our model assumes that the particle energy gain is due 
to both systematic and stochastic particle acceleration, 
while the particle energy losses are due to synchrotron 
emission and adiabatic expansion. 

Describing the systematic acceleration in terms of elec- 
tric field energy gain, we derive a simple relation for the 
expected particle Lorentz factor, 7 ^ lO'^-lO^. Thus, for 
a typical GRB magnetic field of ^10** G, a plasma den- 
sity of ~ 10^ cm^'^,and a beaming factor ~ 50 - 100 the 
expected synchrotron peak energy E'^ is ~ 0.3 MeV as 
found for the observed distribution (e.g. Kaneko et al. 
2006). This may explain the non uniform time resolved 
Ep distribution of the GRB SED peaking around a char- 
acteristic value. 



Following the assumption that systematic and stochas- 
tic acceleration mechanisms balance the synchrotron ra- 
diative cooling, and the adiabatic expansion loss is the 
main process governing GRB spectral evolution during 
the decay phase of individual pulses, we derive that the 
expected intrinsic scaling relation between the height of 

3/2 

the SED Sp and its peak energy Ep is Sp oc Ep' , which 
agrees with the observed HIC for single pulses in long 
GRBs. 

Finally, on the basis of our assumptions, we note that 
the adiabatic losses do not change the shape of the SED 
during the prompt emission. We showed that this is con- 
sistent with the spectral behavior of the decay phase of 
single pulses during long GRB prompt emission as in 
the cases of GRB 941026 and 950818, for which we did 
not detect any large variation in the spectral curvature 
throughout the spectral evolution. 

F. Massaro is grateful to R. Preece for comments 
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Fig. 3. — The time resolved best fit models of the decay phase for 
GRB 921123 performed with the log-paraboIic function. The time 
interval between each SED is 0.5 s, going from the brigthest yellow 
to the faintest red. The cyan dashed line is the time integrated 
best-fit model. 
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